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Abstract 



We show that there are two different families of (weakly) orthog- 
onal polynomials associated to the quasi-exactly solvable Razavy po- 
tential V(x) = (C cosh2x-M) 2 (C > 0, M € N). One of these families 
encompasses the four sets of orthogonal polynomials recently found by 
Khare and Mandal, while the other one is new. These results are ex- 
tended to the related periodic potential U(x) = —(( cos2x — M) 2 , 
for which we also construct two different families of weakly orthog- 
onal polynomials. We prove that either of these two families yields 
the ground state (when M is odd) and the lowest lying gaps in the 
energy spectrum of the latter periodic potential up to and including 
the (M — l) th gap and having the same parity as M — 1. More- 
over, we show that the algebraic eigenfunctions obtained in this way 
are the well-known finite solutions of the Whittaker-Hill (or Hill's 
three-term) periodic differential equation. Thus, the foregoing results 
provide a Lie- algebraic justification of the fact that the Whittaker- 
Hill equation (unlike, for instance, Mathieu's equation) admits finite 
solutions. 
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1 Introduction 



The one- dimensional quantum mechanical potential 

V(x) = (C cosh2x-M) 2 , (1.1) 

where ( and M are positive real parameters, was first studied by Razavy |I| . 
For M > (, the above potential (to which we shall henceforth refer as the 
Razavy potential) is a symmetric double well. This type of potentials has 
been extensively used in the quantum theory of molecules as an approximate 
description of the motion of a particle under two centers of force. In particu- 
lar, the Razavy potential has been proposed by several authors as a realistic 
model of a proton in a hydrogen bond, The potential (|1.1|) has also 

been recently used by Ulyanov and Zaslavskii, ||, as an effective potential 
for a uniaxial paramagnet. 

Razavy showed that when M is a positive integer the lowest M energy 
levels of the potential ( |1 . 1|) (with their corresponding eigenfunctions) can be 
exactly computed in closed form. The Razavy potential is thus an example of 
a quasi-exactly solvable (QES) potential, for which part (but not necessar- 
ily all) of the spectrum can be computed exactly. A very important class of 
QES potentials, that we shall call algebraic in what follows[|, are character- 
ized by the fact that the corresponding quantum Hamiltonian is an element 
of the enveloping algebra of a finite-dimensional Lie algebra of differential 
operators (the so-called hidden symmetry algebra) admitting a finite- 
dimensional invariant module of smooth functions. That such a potential is 
QES follows immediately from the fact that the finite-dimensional module of 
the hidden symmetry algebra is obviously left invariant by its enveloping al- 
gebra, and in particular by the Hamiltonian. Therefore, a number of energy 
eigenvalues and eigenfunctions equal to the dimension of the invariant mod- 
ule can be computed algebraically, by diagonalizing the finite-dimensional 
matrix of the restriction of the Hamiltonian to the module. 

One-dimensional algebraic QES potentials were studied as such for the 
first time by Turbiner || , who used as hidden symmetry algebra a realization 
of st(2, R) in terms of first-order differential operators. These potentials 
were then completely classified by Gonzalez-Lopez, Kamran and Olver @,10]. 



lr There is, unfortunately, no clear consensus in the literature regarding this terminology. 
The term "quasi-exactly solvable potential" is, we believe, originally due to Turbiner and 
Ushveridze, 0, who used it to refer to what we have just called algebraic QES potentials. 
However, in the last couple of years there has been a growing tendency to use the adjective 
"quasi-exactly solvable" for any potential, be it algebraic or not, part of whose spectrum 
can be exactly computed. We have preferred in this paper to adhere to this increasingly 
more common usage to avoid confusion. 
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There are exactly ten families of one-dimensional algebraic QES potentials, 
five of which are periodic and the remaining five all have point spectrum. In 
all cases, the hidden symmetry algebra is again s[(2,R). 

Recently Bender and Dunne |^TJ associated a family of (weakly) orthog- 
onal polynomials to the class of algebraic QES potentials given by 



V{x) 



(4s-l)(4s-3) 
Ax 2 



2(2s + 2J -l)x 2 + x e ; se 



J e N. 



This construction was immediately extended by the authors of this paper 
to virtually all one-dimensional algebraic QES potentials in ||12|| . Krajew- 



ska, Ushveridze and Walczak |T3| proved that a set of weakly orthogonal 
polynomials can be constructed explicitly for any (not necessarily algebraic) 
QES Hamiltonian tridiagonalizable in a known basis. Khare and Mandal 
have constructed two families of non-orthogonal polynomials associated to 
a pair of non-algebraic QES potentials, ]14|]Q- It is important to note that 
the family of polynomials associated to a given Hamiltonian is not unique, 
but depends on the type of formal expansion defining the polynomials. It 
is therefore conceivable that one could obtain orthogonal polynomials in the 
examples studied in [14j] by considering different expansions. 

The Razavy potential has been recently revisited by Khare and Man- 
dal, and Konwent et al., The former authors, who were mainly 
interested in the properties of the associated polynomial system, introduced 
four different sets of polynomials {Pk{E)}^L and {Ql{E)}^L (e = 0, 1) for 
the Razavy potential ( p..l|) through the formulae 



i/>e(x) 



-D' /2 5: 



fc=0 



PRE) fz + i 
(2k)\ 



and 



if> E (x) 



k=0 



Ql(E) 
(2k + l)\ 



z + 1 



where 



z = cosh 2x, (1.2) 

2 Some of the formulae for the polynomials associated to these potentials contain errata. 
Indeed, the change of variable (8) should read t = (y + e 2 ) 1 / 2 , and the factor multiplying n 
in the coefficient of Q n (s) in the recursion relation (11) should be e 2 (4s + 1), this affecting 
formulae (12) and (13). Likewise, in formula (21) the term (a + b + c — |) 2 should be 
{a + b + c + n- §) 2 . 
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and ipE denotes a formal (i.e., not necessarily square- integrable) eigenfunction 
of 

H = -d 2 x + (C cosh 2x - Mf 

with eigenvalue E and parity^ (— Without loss of generality, we shall 
choose the usual normalization 

P*(E) = QUE) = 1. 

Imposing that (H — E) 4>e{x) — one easily shows that each of the four 
sets {Pk(E)}^L and {Q e k (E)}^L (e = 0, 1) satisfies a three-term recurrence 
relation of the appropriate form (see ( |2.12| ) and Ref. |T7[), and therefore 
forms an orthogonal polynomial system with respect to a suitable Stieltjes 
measure. 

We thus have four seemingly unrelated sets of orthogonal polynomials 
associated to the Razavy potential ( |1 . 1| ) . This is surprising, since in all the 
previous examples only one set of orthogonal polynomials was constructed 
for each QES potential considered. One of the objectives of this paper is 
precisely to explain how these four sets of orthogonal polynomials arise. The 
key to this explanation is the fact (not taken into account in [15|]) that the 



Razavy potential is not just QES, but algebraic. More precisely, we shall show 
in Section ^| that the Razavy potential can be written as a polynomial in the 
generators of a suitable realization of s[(2,M) in two different ways. Using 
the constructive method explained in ||12|| , these two different realizations of 
the Razavy potential as an algebraic QES potential give rise to two different 
families of orthogonal polynomials. One of these two families encompasses in 
a natural way the four sets of orthogonal polynomials of Khare and Mandal's. 
In fact, all the properties of these four sets found in |TJ|] (weak orthogonality, 
factorization, etc.) are immediate consequences of the general properties of 
the system of orthogonal polynomials associated to an algebraic QES poten- 
tial developed in our previous paper jE|. The second realization of ( |1.1| ) as 



an algebraic QES potential yields yet another set of orthogonal polynomials 
different from the four sets found by Khare and Mandal. The properties of 
this family, which again follow from the general theory developed in |T2[ , are 



in many respects simpler than those of the first family. For example, the 
moment functional associated to the second family is positive semidefinite, 
while this is not the case for the first family. All of these facts make, in our 
opinion, the second family more convenient in practice for finding the exactly 
computable energy levels of the Razavy potential. 

3 Where (z - l) 1 / 2 should of course be interpreted as \72sinha;. 
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In Section |3] we study the trigonometric version of the Razavy potential, 
given by 



U(x) = — (C cos 2a; — M) 2 . (1.3) 

This potential, which is a simple model for a one-dimensional periodic lattice, 
appears in Turbiner's list of QES one-dimensional potentials, ||, and was also 
touched upon by Shifman, ]18| (in the particular case in which M is an odd 
positive integer). Ulyanov and Zaslavskii, ||, have related the trigonometric 
Razavy potential ( |1.3| ) to a quantum spin system. The potential Q1.3Q has 
also recently appeared as the coupling term between the inflaton field and 
matter scalar fields in theories of cosmological reheating after inflation with 
a displaced harmonic inflaton potential, [|19f| . 

The trigonometric Razavy potential (|1.3|) is the image of the hyperbolic 



Razavy potential ( |1.1| ) under the anti-isospectral transformation x i— > ix 
E i — > — E, recently considered by Krajewska, Ushveridze and Walczak |20| 



It is therefore to be expected that the properties of the polynomials associ- 
ated to this potential are analogous to the corresponding properties for the 
hyperbolic Razavy potential ( |1 . 1| ) . That this is indeed the case is shown in 
Section |3], where we prove that the potential (|1.3| ) can be realized in two 
different ways as an algebraic QES potential. As in the hyperbolic case, each 
of these two different realizations gives rise to a family of orthogonal polyno- 
mials. For each positive integer value of M it is possible to exactly compute 
M eigenf unctions (with their corresponding energies) of the trigonometric 
Razavy potential by purely algebraic procedures. It was previously known 
on general grounds that the energies of these M eigenf unctions must be 
boundary points of allowed bands (or, equivalently, gaps) in the energy spec- 
trum of the periodic potential ( |1.3j ). In Section |3] we investigate the exact 



position of these boundary points in the spectrum of (.L3). We will show 
that, if the gaps in the energy spectrum are numbered consecutively in order 
of increasing energy, these points yield precisely the ground state (when M 
is odd) and the lowest [4^-1 gaps^of the same parity as M — 1. For instance, 
if M = 4 we obtain the first (lowest) and the third gap, whereas for M = 5 
we get the ground state and the second and fourth lowest gaps (see Fig. |l]). 

The paper ends with a discussion of the above results in the context of 
the classical theory of Hill's equation. We show that the algebraic eigenfunc- 
tions constructed in this paper are precisely the so called finite solutions of 
the Whittaker-Hill (or Hill's three-term) equation. In fact, our analysis pro- 
vides a Lie-algebraic explanation of why the Whittaker-Hill equation admits 
finite solutions at all. Indeed, from our point of view this is just a simple 



4 We denote by [x] the integer part of the real number x. 
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M 



M 



Figure 1: Structure of the energy spectrum of the potential ( |1.3| ) 
for M = 4 and M = 5. The thick horizontal segments repre- 
sent the allowed energy bands. The solid circles stand for the 
algebraically computable energies, which determine the gaps rep- 
resented by thin solid lines. The remaining boundaries of the 
allowed bands, and hence the energy gaps shown as dotted lines, 
cannot be exactly computed. 



consequence of the fact that the Schrodinger operator with potential is 
algebraically QES. 

2 The hyperbolic Razavy potential 

We shall show in this section that the hyperbolic Razavy potential can 
be expressed in two different ways as an algebraic QES potential. From 
these two representations we shall derive two different families of associated 
orthogonal polynomials, whose properties we shall discuss. 

Consider, in the first place, the second non-periodic algebraic QES po- 
tential listed in Ref. [K| (p. 127), given by 

V(x) = A cosh 2 \Jv x + B cosh \fvx + C coth \pu x csch \Jv x + D csch 2 \fvx . 

(2.1) 

where the coefficients A, B, C, D can be expressed in terms of four parameters 
b, c, d G E and n e N U {0} as follows (see 0, Eq. (5.11)): 

A=^, B = ±[c + (n + l)v), C= b -±^[c-(n+l)u] 



D = — 
Av 



+ df+ (c - (n + l)v) 2 - v 2 . (2.2 



The hyperbolic Razavy potential is of the form Q2.1D (up to an inessential 
additive constant) provided that 

v = 4, C = D = 0, A = ( 2 , B = -2MC. 
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d 


M 


a 


V 


4C 


2n + 1 


1 





4C 


2n + 3 


-1 





4(C-1) 


2n + 2 





1 


4(C + 1) 


2n + 2 





-1 



Table 1: Values of d and M corresponding to the four solutions 
of Eqs. (pl)-(pl). 



Using Eq. ( j2.2| ) we obtain the following system in the parameters b, c, d 
and n: 

b 2 = lQ( 2 , (2.3) 
6 [c + 4(n+l)] = -16MC, (2.4) 
(6 + d) [c - 4(n + 1)] = 0, (2.5) 
(6 + rf) 2 + [c-4(n+ l)] 2 = 16. (2.6) 

From Eqs. ( |2.3| ) and (|2.4[ ) and the normalizability condition b < (see 0, 
Eq. (5.14)), we get 

S = -4C, c = 4(M-n-l). 

Substituting these into the remaining conditions ( p.5| ) and ( |2.6| ) we are led 
to four different solutions for d and M, which may be written in a unified 

way as 

M = 2(n+ 1) -a, c=4(n + l-a), d = A{(-r]), 

where the parameters a and rj are given in Table II]. 

It follows from the general discussion in Ref. || that the change of vari- 
able (|1.2|) and the gauge transformation determined by 

fa) = (z- l^-'-^iz + i)i(i-^) e -^/2 (2 _ 7) 
map the Razavy Hamiltonian into an operator H g (the gauge Hamiltonian) 
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belonging to the enveloping algebra of the realization of sl(2, M) spanned byQ 



ft 

J_ = d z , J = zd z - - , J + = z 2 d z -nz. {21 



Indeed, 



H g {z) = — ■ [-d 2 x + (Ccosh2x - M) 2 ] • /2( 



where 



x= i arccosh z (2.9) 

-4 (j 2 - J 2 - CJ+ + (n + 1 - a) Jo + (C " rf)J-) - c* , 

c, = -3(n + l) 2 + 2(n + l)a + 2r/C - C 2 - 

According to the general prescription of Ref. [E2| , the formal solutions of the 
gauged equation 

(H s -E)xe = 0, (2.10) 

are generating functions for a set of orthogonal polynomials. More explicitly, 
inserting the expansion 

into ( |2.10| ), we readily find that the coefficients Pk{E) satisfy the three-term 
recursion relation 

P k+1 {E) — (E — 6 fc ) P k (E) - a fc Pk-xiE), k>0, (2.12) 

where 

a fc = 16(k(2k - a + rj)(k - n - 1) , 

6 fc = -4A;(A; + 1 - a + 2() + (2n + l)(2(n - a) + 3) + C(C - 2?7 + 4n) . 

(2.13) 

If we impose the condition Pq(E) = 1, the coefficients {Pk^E)}^^ form a 
set of weakly orthogonal (monic) polynomials. Therefore, we can construct 
two sets of weakly orthogonal polynomials for each value of M by choosing 
suitable values for a, r\ and n (for M = 1 there is only one set). These 

5 The operators J a (a — ±,0) (and any polynomial thereof) preserve the space V n 
of polynomials in z of degree at most n. Moreover, any fc-th order differential operator 
(fc < n) preserving P n may be expressed as a fc-th degree polynomial in the generators 

Ja, 1311. 
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two sets coincide exactly with the sets {P k 1 ~ e (E)} < j? =0 and }fcL (with 

e = for M even and e = 1 for M odd) studied by Khare and Mandal 
15 ] . From now on, we shall use when necessary the more precise notation 



in 



to denote the orthogonal polynomials defined by Eqs. (|2.12| )-( p.l3| ). For 
instance, if M = 4 we have n = 1, a — and rj = ±1 (see Table |l|). When 
rj = 1 the first three polynomials are 

P° + (E) = 1, 

P? + (E) = E-( 2 -2(- 15, (2.14) 
P$ + (E) = E 2 - 2(C 2 - 2C + 11)£ + C 4 - 4C 3 + IOC 2 - 28C + 105, 
whereas for i] = — 1 we obtain 

A 0- (£) = E - C 2 - 6C - 15, (2. 15) 

P 2 °-(£) = E 2 — 2(C 2 + 2C + 11)£ + C 4 + 4C 3 + IOC 2 + 28C + 105. 

The polynomials (|2.14| ) and (|2.15|) reduce, respectively, to the polynomials 
Q° k (E) and P£(E) (k = 0, 1, 2) in formulae (2.18) and (2.17) of Ref. p| . 



The coefficient a k given by ( |2.13| ) vanishes for k = n + 1, and therefore 
Pk(E) with k > n+ 1 factorize as P n+ i + j(E) = Qj(E)P n+ i(E), j > 0, where 
{Q^JLq also form a set of (monic) orthogonal polynomials. If Ej is a root 

of the polynomial P n+ i(E) the series ( |2.11| ) truncates at k — n, and thus 
Ej belongs to the point spectrum of the Razavy Hamiltonian. For example, 
if M = 4, the roots E , E 2 of p2 + (E) are the energies of the ground state 
and the second excited state of the Razavy potential, while the roots E%, E 3 
of p2~(E) correspond to the first and third excited levels. The rest of the 
spectrum cannot be computed algebraically. 

The other usual properties which characterize weak orthogonality - 
vanishing norms, finite support of the Stieltjes measure associated to the 
polynomials, etc., fTT|-|T3| — are also satisfied by the polynomials {Pfc(-E)}£i . 
In particular, if E^ (k — 0, . . . , n) are the (different) roots of P n+ i(E), the 
moment functional associated to the polynomials is 

n 

C = J2 UJ kS(E-E k ), (2.16) 
where the coefficients u k are determined by 

n 

^2 Pi(E k )u k = 5i , l = 0,...,n. 

k=0 
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It was observed in []ToJ that not all the coefficients oj k corresponding to the 
polynomials P| and Q\ are positive. This is in fact a direct consequence of 
the following general property of an orthogonal polynomial system satisfying 



a recursion relation of the form ( |2.12| ) 



Proposition 2.1 The coefficients u k of the moment functional (|2.16|) are 

positive for all k — 0, . . . ,n if and only if a k > for < h < n and b k is real 
for < k < n. 

Proof. The "if" part was proved in [00. Let u k > for k — 0, . . . , n. Then 
C(p 2 ) > for any non- vanishing real polynomial of degree at most n. It 
follows that C is positive definite in P 2 n = {P £ C[E] : degP < 2n}, for if 
P ^ P2n is a nonzero real polynomial which is non-negative for all £gI, 
then P = p 2 + q 2 for real polynomials p and q in V n , and thus C(P) > 0. Since 
£ is positive definite in P 2n , the moments /i^ = C(E k ) with k = 0, . . . , 2n 
are positive for even k and real for odd k, fT7|| . Multiplying the recursion 
relation ( |2.12j) by P fc and applying £ we find that 

C{EPl) - b k C(P 2 k ) = . (2.17) 

Taking /c = 0, we conclude that bo = fii/fio is real. Therefore Pi = E — b Q 
is real and a\ = C(P 2 ) > 0. By induction, if bj-i G M. and a, > for 
j = 1, . . . , k < n, then P k is real, and from ( p,17|) we deduce that b k G M. 
Then P k +i is real, and 

fc+i 

o<£(p fc 2 +1 )=n^ 

implies that a k+ i > 0. Q.E.D. 

Note that the coefficients a k given by ( |2.13| ) are negative for 1 < k < n 
and therefore u k cannot be positive for all k = 0, . . . , n. 

Consider, in the second place, the third non-periodic algebraic QES po- 
tential given in Ref. [[H]] (p. 127), namely 

V(x) = Ae 2 ^ x + Be VVx + Ce~^ x + De- 2V " X , (2.18) 

where the coefficients A,B,C,D can again be expressed in terms of four 
parameters b,c,d G M and n G N U {0} as (see 0, Eq. (5.11)): 

A=^ P = A[c+(n + lH, C = ±[c-(n+l)u], P = ^ 

(2.19) 



3 Note that the proof only requires b k to be real for < k < ri- 



ll 



The potential ( |2.18| ) reduces to the hyperbolic Razavy potential ( |1.1| ) (up to 
an additive constant) provided that 

c 

z, = 4, A = D = -j, B = C = -(M. 

Taking into account the normalizability conditions b < and d > 0, (see |§, 
Eqs. (5.24) and (5.25)), we get the unique solution 

b = -2C, c = 0, d = 2(, M = n + 1. 

In this case, the change of variable 

z = e 2x (2.20) 

and the gauge transformation generated by 

= z^F-e-te+k) (2.21) 

map the Razavy Hamiltonian into a differential operator H g quadratic in the 
generators (|2.8|) , namely 

H g (z) = -4J 2 + 2C J+ - 2C J- - £* , (2.22) 

where 

c* = -(n + l) 2 -C 2 - 



Following the general treatment of ||12|| , if we insert the expansion 

fe=0 v s; 

into the spectral equation for if g , the coefficients Pk(E) are easily found to 
satisfy a three-term recursion relation of the form fl2.12|) , with coefficients 



a k = Ak(n + 1 - k)( 2 , 
b k = 4k{n - k) + 2n + 1 + ( 2 . 



(2.24) 



Taking Pq{E) = 1, we obtain yet another family of weakly orthogonal (monic) 
polynomials {Pk(E)}^L associated to the Razavy potential fll.lj) . For in- 
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stance, if M — 4 the first five polynomials are 

Po(E) = 1, 
P l (E)=E-( 2 -7, 

P 2 {E) = E 2 — 2(C 2 + 11)E + C 4 + IOC 2 + 105, 

P 3 (E) = E 3 - (3C 2 + 37)E 2 + (3C 4 + 46( 2 + 435)£ - C 6 - 9( 4 - 143C 2 - 1575, 
P 4 (E) = E 4 — 4(C 2 + 11)£ 3 + 2(3C 4 + 46C 2 + M7)E 2 - 4(C 6 + 13( 4 + 159( 2 
+ 1155) E + C 8 + 4C 6 + 86C 4 + 1316C 2 + 11025. 

(2.25) 

Note that the polynomial Pi(E) is the product of the polynomials p2 + (E) 
and p2~(E) given in (|2.14|) and (|2.15|) . Therefore, the algebraic levels can 
be also obtained as the the roots E , . . . ,E 3 of P±{E). 

In general, if M is even the polynomial Pm{E) factorizes into the product 
of the polynomials P^j 2 (E) associated to 77 = ±1 (see Table |l|). Alterna- 
tively, if M is odd, Pm(E) factorizes into the product of the polynomials 
(m±i)/2 associated to a = ±1. The algebraic energy levels of the Razavy 
potential ( |1.1|) can thus be computed in a unified way as the roots of Pm- 
On the other hand, the algebraic eigenf unctions can be written as 

ip E (x) = h(z)xe(z)\ z=z(x) , 

where fi(z), Xe(z) and the change of variable z = z(x) are given by ei- 
ther (|T2|), Q and (FTTH) , or (|2T20|) , ( ^21| ) and (|2T23|) . 



The polynomials {Pk}^L verify the usual properties associated to their 
weak orthogonality. However, unlike the previous family {Pk}kLo, the coef- 
ficients hk of the recursion relation are positive for < k < n. It follows 
from Proposition |2.1| that the coefficients u>k of the corresponding moment 
functional £ are positive for all k — 0, . . . , n, i.e., £ is positive semidefinite. 

Before finishing this section, let us emphasize that the Razavy Hamil- 
tonian admits two different gauged forms H g and H g inequivalent under 
the action of the projective group on the enveloping algebra of the genera- 
tors ( p.8|) , . This does not contradict the fact that H g and H g are equivalent 



under a change of variable and a gauge transformation (since they are both 
equivalent to the Razavy Hamiltonian) . Indeed, the transformation relating 
H g and H g , 

HJi) = I ^ H.(z) < ^ 



is certainly not projective. 
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3 The trigonometric Razavy potential 

3.1 The orthogonal polynomial families 

We shall study in this section the trigonometric Razavy potential 

U(x) = -{( cos 2x - Mf , (3.i; 



which can be obtained from the hyperbolic Razavy potential ( |1 . 1|) applying 
the anti-isospectral transformation x > i x, E 1— > —E. In other words, if)(x) 
is a solution of the differential equation 

[-91 + V(x)] 7p(x) = Eip(x) (3.2) 

if and only if 

</>(x) = ij){ix) (3.3) 

is a solution of 

[-dl + U(x)] </>(x) = -E<l>(x). (3.4) 

Just as in the hyperbolic case, we see by inspection that the trigonometric 
Razavy potential can be expressed as an algebraic QES potential in two 
different ways. Indeed, the potential ( |3.1| ) is a particular case of two entries 
in the table of periodic one-dimensional QES potentials given in Ref. flTOfl : 

case 4, 

U(x) = A sin 2 \fvx + B sin \fvx + C tan y/ux sec \fvx + D sec 2 \fux 

(3.5) 

for 

z/ = 4, A = -(\ B = 2M( } C = D = (3.6) 

(after performing the translation x 1— > x — 7r/4; notice that one-dimensional 
QES potentials where classified in Ref. [|K| only up to an arbitrary transla- 
tion), and case 5, 

U(x) = A cosAy/ux + B cosZyfvx + C sva.2y/i}x + D smAy/vx (3.7) 
with 

c 2 

v = l, A=-y, B = 2M(, C = D = 0. (3.8) 
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The four parameters A, B, C, D appearing in the table of periodic QES po- 



tentials of Ref. [10] are not independent, but must satisfy a single algebraic 
constraint; it can be verified that the choices of parameters ( |3.6| ) and ( |3.&| ) do 
indeed satisfy this constraint. Let us also note at this point that, although 
the representation (|3~5| )-(|3~6|) was known to Turbiner, ||, and Shifman, |18| . 
the second representation appears to be new. 

Let us now construct the systems of weakly orthogonal polynomials as- 
sociated to the representations ( |3.5| ) and ( |3.7| ) of the trigonometric Razavy 
potential QES potential. 

Consider, in the first place, the representation Q3.5p . Instead of proceeding 
directly, along the lines sketched in the previous section for the hyperbolic 
case, we shall exploit the fact that the trigonometric and the hyperbolic 
Razavy potentials are related by the anti-isospectral transformation 
We saw in the previous section (Eq. (|2.9[)) that 



-d 2 x + V(x) = ft(z) ■ H g {z) 



1 



fi(z) 



=cosh 2x 



performing the change of independent variable x ^ ix we obtain 

1 



or, equivalent ly, 



di + V{ix) = fl(z) ■ H g (z) 



-d 2 x + U(x) = - fi(z) ■ H g (z) 



jl{z) 



z=cos 2x 



fi{z) 



2=cos 2x 



Thus, the gauge Hamiltonian if| er associated to the potential U(x) in this 
case is simply 



m e *(z) 



-HJz) 



(3.9) 



Since (cf. Ref. [0, Eq. [41]) the coefficients and bk defining through 
Eq. (p.!2[) the orthogonal polynomial system associated to a QES potential 
are, respectively, quadratic and linear in the gauge Hamiltonian, it follows 
that the recurrence relation satisfied by the orthogonal polynomial system 



Pl er (E) [> associated to the representation (|3.5|)-(| 



is 



P^\(E) = {E + b k ) Pr(E) - a k Pp-(E), k > 
Comparing with ( p,12|) we immediately obtain the relation 

k > 0. 



pr(E) = (-iyPk(-E) 



(3.10) 
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It can be easily verified through a routine calculation similar to the one 
performed in the previous section that the general procedure described in 
Ref. |T2[ to construct the orthogonal polynomial system associated to a QES 
potential, when applied to the representation ( |3.5| ) of the potential ( |3.1| ), 
does indeed yield the result (|3.10|) . 

Let us now turn to the representation ( [3.7|) . We cannot directly apply 
the previous reasoning in this case, since the composition of the change of 
coordinate z = e 2x with the anti-isospectral mapping x i— > ix leads to the 
complex change of variable z = e 2lx , while the correct one for this case is 

(cf. H) 

z = tanx . (3-11) 



It is therefore easier to apply, as in the previous section, the general procedure 
described in [0. Using the techniques explained in Ref. |J, we readily find 
the following expression for the coefficients A,B,C,D in (|3.7|) in terms of 
four independent parameters b,c,d GK and nGNU {0}: 



A 



C 



1 



1 



32 v 
1 



B = — 

8u 



16 v 



Comparing 



(b-c- d){b + c-d), 

OA 

c(b + d) +Au(b-d)(n+ 1)1 , D 

with ( p. 12 ) we readily obtain 
b = d = 0, c = 4(, n = M - 1 . 



d 2 - b 2 + Avc(n + 1) 



(d-b). (3.12) 



(3.13) 



From ( |3.13|) it follows (cf. Ref. 0) that the change of variable ( |3.11|) and the 
gauge transformation determined by 



/P er (z) 



^ 2 + l)- 



2 e 



(3.14) 



map the trigonometric Razavy Hamiltonian into the gauge Hamiltonian 



#r r w = - 

in the sense that 



J\ + 2 J 2 + J 2 + 4C J + -(M 2 + 2C 2 + 1) 



-d 2 x + u(x) = ~^{z) ■ m<*{z) 



flP ei {z) 



(3.15) 



(3.16) 



As in the previous section, the differential operators J e (e = 0, ±) appearing 
in (|3.15|) are defined by (|2.8|), with n — M — 1. The orthogonal polynomial 
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system < P£ eT (E) > for this case is generated by expanding an arbitrary 



fc=0 



solution x^ er (-2) of the gauged equation 

E 



H? er (z) 



X V e(z) 



(3.17) 



in the formal power series (cf. Ref. |12[, Eqs. (29) and (40)) 



xT(*) 



[z + i 



00 ( — A\ k 



k=0 



(2() k k\ 



z — I 



z + i 



(3.18) 



From (|3~T6D and ( gl7|) it follows that 

is a formal solution of the Schrodinger equation 

[-d 2 x + U(x)} fc(x) = Ej%'(x) 



z=tan x 



(3.19) 



(3.20) 



Applying the change of variables x \— > —ix we deduce that ip^(—ix) satisfies 
the dual equation 



l-d 2 x + V(x)} r_%(-ix) = Eft%{-ix) . 
On the other hand, a direct calculation shows that 



'0 



per . 



-IX) 



e 2 



k=0 



(2() k k\ 



(3.21) 



(3.22) 



where jl(z) is defined by ( |2.21[ ). From ( |3.21| ) and ( |3.22| ) it follows that the 
formal power series 



k=0 



(2C) fc A;! 



(3.23) 



must be proportional to the function Xe{z) generating the orthogonal poly- 
nomial system associated to the the representation ( |2.18|) - ([2.19|) of the hy- 
perbolic Razavy potential as an algebraic QES potential (see Eq. (|2.23|) ). 
Since both power series ( |3.23| ) and ( |2.23| ) have constant term equal to one, 
they must be equal. Equating their coefficients we obtain the relation 



Pr(E) = (-l) h P k (-E), k>0. 



(3.24) 
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3.2 The band spectrum 

The results of the previous subsection imply that, as in the hyperbolic case, 
the trigonometric Razavy potential is algebraically QES when MeN (see, 
for instance, Eq. ( ft.l3| )). We shall see in this subsection how this fact can 
be used to exactly compute a certain number of gaps in the band energy 
spectrum of the potential ( |3.1| ), and shall furthermore determine the location 
of these gaps in the spectrum. 

We start by briefly recalling certain well-known facts about the energy 
spectrum of a Schrodinger operator 

H = -dl + U(x) , (3.25) 

where U is a continuous periodic function of period a > 0: 

U(x) = U(x + a), Vxel. (3.26) 

A real number E belongs to the spectrum of H if the differential equation 

(H-E)ip = (3.27) 

has a bounded nonzero solution tp(x). It can be shown, P3|, |2"3| , that the 
spectrum of the operator ( |3.25| )-( |3.26| ) is the union of an infinite number of 
closed intervals (energy bands) [Eq, Ei], [E 2 , Ei], [E 2 , E 3 ], . . . , where 

E < Ex < E 2 < E l < E 2 < E 3 < E A < . . . 

and 

lim Ek = lim Ek = +oo . 

The gaps in the energy spectrum are thus the (possibly empty) open intervals 

(Ei, E 2 ) , (Ei, E 2 ) , . . . , (Ek, E k+ i) , (E k , E k+ i) , . . . ; k — 1,2, ... . 

The numbers Ek (k = 0,1,2,...) are characterized by the existence of a 
nonzero a-periodic solution ipk of the differential equation (|3.27 ) for E = Ek- 
The latter condition is clearly equivalent to ipk being a nonzero solution of 
the Sturm-Liouville problem with periodic boundary conditions 

(H - E)ip(x) = 0, 0<x<a; 

= V'(0) = ^'(o), (3.28) 

with eigenvalue E = E k . Note that Ek = Ek+i if and only if ( |3.28|) has two 
linearly independent solutions. Likewise, the numbers E k (k = 1, 2, . . . ) are 
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characterized by the fact that the differential equation ( |3.27| ) with E = E k 
possesses a nonzero anti-periodic (2a-periodic) solution ip k , that is, a solution 
ip k such that 

ip k (x + a) = —tp k (x) , el. 
Equivalently, ipk is a solution of the Sturm-Liouville problem 

(H - E) if>(x) = , 0<x<a; 

^,(0) = -ip(a) , V'(0) = -^'(a) (3.29) 



with eigenvalue E = E k . Finally, it is shown in Ref. [2S| that for = 
0,1,2,... the eigenfunction ip k has exactly k + 7r(/c) zeros in the interval 
[0, a), where 

*(*) = i [i + (-i) fc+1 ] 

is the parity of k. A straightforward adaptation of Ince's proof shows that 
ijjk has exactly k + 7r(k) — 1 zeros in [0, a), where now k = 1, 2, . . . . 

Let us now turn to the trigonometric Razavy potential ( |3.1| ), which we 
know from the previous discussion to be algebraically QES for MeN. For 
convenience, we shall use in what follows the representation ( |3.7| )-( pT8| ) of the 
trigonometric Razavy potential as an algebraic QES potential. From (|3.24j) 
it follows that the polynomials P^ er (E) satisfy the recurrence relation 

^ r i(£) = (E + h) Pr(E) - a k P£\(E) , k > , 

where the coefficients a k and b k are given by Q2.24Q with n = M — 1. In 
particular, the coefficient au vanishes identically, which in turn implies that 
if 6j is a root of the polynomial P^ r (E) then the series ( |3.18| ) truncates at 
k — M — 1, and therefore the function if)f° r (x) given by ( |3.19| ) is a regular, 
bounded solution of the Schrodinger equation ( |3.20| ). Since it can be shown 



12| that the polynomial PYT has exactly M different real roots, we can 



M 

algebraically compute M solutions of the Schrodinger equation ( |3.20| ) of the 
form (cf. Eqs. frf| ), ( gig ), and (|3l9|) ) 



^(x) = e* ™*2*-W-T)x (f . ( e 2ix^ ^ 1 <j < M , (3.30) 

where ej is any of the roots of P^ r (E) and tpj is a polynomial of degree at 
most M — 1. From the latter equation we find that 

^(x + n) = (-l) M -^(x) } (3.31) 

and therefore the M algebraically computable solutions ( |3.30| ) are 7r-periodic 
for M odd, and anti-periodic for M even. Furthermore, from ( |3.30|) it also fol- 



lows that each of the M algebraic eigenfunctions ip^ ev has at most M— 1 roots 
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in the interval [0, 7r). Thus for M even the algebraic eigenf unctions ( |3.30| ) co- 
incide with the M lowest anti-periodic eigenfunctions ^k{x) (k — 1, . . . , M), 
and the algebraically computable energies — the zeros ej of the critical poly- 
nomial P^ r {E) - - are the lowest M "anti-periodic" energies Ei, . . . , E M . 
Likewise, for M odd the algebraic eigenfunctions are the lowest M periodic 
eigenfunctions ipk{x) (k = 1, . . . , M), and the corresponding algebraic ener- 
gies are the M lowest "periodic" energies E , . . . ,Em-i- Since the gaps in 
the energy spectrum of a periodic Hamiltonian are limited by two energies of 
the same type (periodic or anti-periodic), this means that the knowledge of 
the algebraic energies allows us to exactly compute a certain number of gaps 
in the spectrum. More precisely, when M is even then we can algebraically 
compute the first M/2 anti-periodic gaps (E\, E 2 ), (E 3 , £4), . . . , (Em-i, Em) 
in the energy spectrum of the trigonometric the Razavy Hamiltonian or, 
equivalently, the first M/2 odd gaps. Similarly, when M is odd then the 
algebraically computable energies are the ground state E and the first (M — 
l)/2 = [M/2] periodic gaps (E 1} E 2 ), (E 3 , E 4 ), . . . , (E M -2, E M -i) or, what 
is the same, the ground state and the first [M/2] even gaps in the energy 
spectrum. In particular, for M even we can algebraically compute the first 
gap (Ei,E 2 ) in the energy spectrum of the trigonometric Razavy Hamilto- 
nian, while for M odd we can always compute the ground state E . Note 
also that, since the algebraically computable gaps are never consecutive, we 
cannot exactly compute any of the allowed energy bands. Figure |2] shows 
the first five allowed energy bands for the trigonometric Razavy potential as 
a function of the parameter ( for M = 5. 
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Figure 2: Allowed energy bands (in gray) for the potential ( |3.1| ) 
with M = 5 as a function of £. The solid boundary lines corre- 
spond to the algebraically computable (periodic) energies, while 
the gray boundary lines have been obtained numerically 



The differential equation (|3.27f) with the potential ( |3.1|) is well-known in 



the theory of periodic differential equations under the name of Whittaker- 



Hill's equation, or Hill's three-term equation, pq, p7[. It is of interest in 
the latter context mainly because, unlike the much better known Mathieu's 
equation, for certain values of the spectral parameter E it admits so called 
finite solutions, i.e., solutions of the form 

if>(x) =e-2 cos23 V(x), 



where f(x) is a trigonometric polynomial. It follows from Eq. ( |3.30| ) that 



the algebraic eigenf unctions obtained in this section are finite solutions. In 
fact, the converse also holds, namely all finite solutions are algebraic eigen- 
functions. This follows at once from Theorem 7.9 of Ref. which states 



(in our notation) that for each M G N the Whittaker-Hill equation has at 
most [M/2] gaps of periodic (if M is odd) or anti-periodic (if M is even) 
type. For even M, we have shown that there are exactly M anti-periodic 
algebraic eigenfunctions, whose corresponding M eigenvalues define exactly 
M/2 gaps of anti-periodic type. Since, for these values of M, there are also 
exactly M anti-periodic finite solutions and eigenvalues, |26] , it follows from 
the Theorem quoted above that the finite solutions coincide (up to a constant 
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factor) with the algebraic eigenfunctions. A similar argument is valid when 
M is odd. 

In order to compare our findings with the classical theory, it is more con- 
venient to use the representation ( |3.5|) - (|3.6|) . Applying the anti-isospectral 
transformation to Eqs. ( |2.7|) and (|2.11|) of Section H, it follows that the alge- 
braic (unnormalized) eigenfunctions can be classified as follows: 



M even: 



M 



-1 



^( x ) = e -^ cos2x E l cos 2k+1 x, (3.32) 

fc=0 ' /■ 



M 



2 1 P°~(— F) 

^(x) = e-* cos2x sinx \ 7a, cos 2k x; (3.33) 



M odd: 



M-l 

ij,(x) = e -* cos2x k ( n h v cos 2k x, (3.34) 



(2k) 



k=0 

M-3 



^ (2k +1)1 y J 



k=0 

In the above formulae, the polynomials are defined by the recursion 
relation ( [2.12 )-( |2.13| ), and E is one of the algebraically computable ener- 
gies, i.e, —E is a root of the critical polynomials P^j 2 , P^ 2 , Pn^ +1 y 2 ^ 
and P7m-i)/2i respectively. Comparing with the formulae in Section 7.4.1 



of Ref. [26] we easily find that, if ip(x) is an algebraic eigenfunction of one of 
the four types flO^ -flOS!) , then cos2x ip(x) is proportional, respectively, 
to the Ince polynomial Cfffitl, S 2k jt\, Cm-d anc ^ ^m-d where (in the 
notation of Ince, cf. Ref. p6|) E — a^t\, b 2 ^t\, o 2k I _ 1 , or 6^|_ l3 respectively. 

The results of this section can therefore be interpreted as providing a 
deep Lie- algebraic justification for the exceptional fact that the Whittaker- 
Hill equation admits finite solutions. This observation is further corroborated 
by the fact that other periodic Schrodinger equations known to have finite so- 
lutions (in a slightly more general sense) as, for instance, the Lame equation, 
are also algebraically QES, P,|ID|,p5|,p9|. The above results underscore the 



close connection between the existence of finite solutions of Hill's equation, 
on the one hand, and the algebraic QES character of its potential, on the 
other. This remarkable connection certainly deserves further investigation. 
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